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Understanding the role of spin-orbit comling (SOC) has been crucial to controlling magnetic 
anisotropy in magnetic multilayer films [JQ. It has been shown that electronic structure can be 
altered via interface SOC by varying the superlattice structure, resulting in spontaneous magneti¬ 
zation perpendicular or parallel to the plane UM- In lieu of magnetic thin films, we study the 
similarly anisotropic helimagnet Cri/ 3 NbS 2 , where the spin polarization direction, controlled by 
the applied magnetic field, can modify the electronic structure. As a result, the direction of spin 
polarization can modulate the density of states, and in turn affect the in-plane electrical conductiv¬ 
ity. In Cri/ 3 NbS 2 , we found an enhancement of in-plane conductivity when the spin polarization 
is out-of-plane, as compared to in-plane spin polarization. This is consistent with the increase of 
density of states near the Fermi energy at the same spin configuration, found from first principles 
calculations. We also observe unusual field dependence of the Hall signal in the same temperature 
range. This is unlikely to originate from the non-collinear spin texture, but rather further indicates 
strong dependence of electronic structure on spin orientation relative to the plane. 

PACS numbers: 


Despite the fact that its typical energy scale in 3d fer¬ 
romagnetic metals is small compared to other relevant 
scales such as band widths, SOG mixes the nature of 
the spin and orbital components of the Bloch state in a 
nontrivial way and leads to a variety of electrical trans¬ 
port phenomena e.g. the anomalous Hall effect (AHE), 
anisotropic magnetoresistance (AMR), and the planar 
Hall effect. In addition, the recent work in non-collinear 
magnetically ordered states and the related topological 
Hall effect [7H9j not only has renewed the pivotal role of 
SOG through the Dzyaloshinskii-Moriya (DM) interac¬ 
tion (lol - [ll| . but also has presented a possibility to em¬ 
ploy these findings for functional components in magnetic 
devices [3 [11. Non-collinear magnetic ordering is also 
suggested to possibly manifest spin-orbit coining in a 
complex manner, through the DM interaction|12lll5Lll6l|. 
Consequently, the modification of electronic structure by 
spin-orbit coupling is expected to make in-plane electri¬ 
cal transport sensitive to the magnetization orientation 
relative to the plane. 

Cri/ 3 NbS 2 has a layered crystalline structure, in which 
3d transition metal Cr atoms are intercalated in the 
hexagonal 2 iT-type NbS 2 matrix as trivalent ions and 
magnetically order at Tq = 133 K. The ferromagnetic 
layers of Cr^+ lie coplanar with the crystallographic ab- 
planes and the magnetic helix propagates along the c- 
axis with a long pitch of 48 nm, corresponding to 40 unit 
cells da. Its helimagnetic ordering is attributed to the 
DM interaction, which originates from a broken inversion 
symmetry shared by all members of space group P 6322 


With application of a magnetic field {H) along 
the a&-plane, i.e. perpendicular to the helical axis, ferro¬ 
magnetic domains are created between the winds of the 
helix, increasing the length of the magnetic unit cell and 
forming the chiral soliton lattice phase As field is 

increased, all spins become polarized at = 0.18 T. 
Alternatively, when H is applied along the c-axis, i.e. 
along the direction of the helical wave vector, the helices 
smoothly evolve through a conical state. The conical 
angle decreases with increasing H until polarization at 
Hp = 2.5 T (3) [3l. Its electrical conduction is quasi 
two-dimensional: the electrical resistivity measured with 
current flowing along the c-axis is in the order of 10 ^ — 10 ^ 
times larger than with current in the a&-plane 17|, [2l| . In 
these regards, Gri/ 3 NbS 2 displays similar magnetic and 
structural anisotropies as fabricated planar devices. 


So far, SOC in magnetic multilayers or superlattice 
has been studied intensively in terms of engineering the 
electronic structure in order to control the spontaneous 
magnetic anisotropy [1, Q. Here, we demonstrate the 
reverse process in a similarly anisotropic layered system 
Cri/ 3 NbS 2 , such that the spin polarization direction con¬ 
trolled by applied magnetic field alters the electronic 
structure via SOC. For spin polarization along the c- 
axis, this results in enhanced electrical conductivity, by 
increasing the DOS near the Fermi surface, as supported 
by our first principles calculations. This observation is 
the most prominent for temperatures T < Ta', Ta is the 
temperature scale below which spin-disorder scattering 
contributes negligibly to the resistance, and is empiri- 
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cally determined from the temperature dependence of the 
magnetoresistance. Surprisingly, roughly the same tem¬ 
perature scale of Ta was found in the transverse Hall re¬ 
sistivity, and is also consistent with the onset of deviation 
from Bloch’s law for the magnetization [l^. While 
Ta is a crossover temperature and is not sharply defined, 
it corresponds to the energy scale associated with modi¬ 
fications of the electronic structure, and is thus a useful 
quantity to estimate. 

Single crystals of Cri/ 3 NbS 2 were grown by vapor- 
transport (see Ref. 0 for details on growth and sample 
characterization). Crystals from two different batches 
were used in this study, and show slightly different Tcs 
of 120 K and 133 K, but otherwise exhibit very lit¬ 
tle qualitative difference in the T and H dependences 
of both electrical and magnetic measurements. All the 
data presented in this letter are from the sample with 
Tc = 133 K. Samples from both batches were character¬ 
ized with a;-ray diffraction for the crystalline structure 
and small angle neutron scattering to verify the helimag- 
netic ground state [2^. Magnetization was measured in 
a superconducting quantum interference device magne¬ 
tometer. Electrical transport properties were studied us¬ 
ing standard four-probe measurements on samples with 
typical area « 2 x 0.3 mm^ and thickness ~ 20 — 70 
/im. Electrical contacts, made by silver paint, had typ¬ 
ical contact resistance less than 1 H. Rotation of the 
sample in a field was performed by a home-built rotation 
probe, which was inserted in a superconducting split coil 
magnet. For transport measurements, currents of 2 — 4 
mA were applied in the ab plane. The magnetization 
measurements were performed using SQUID. We used 
small wedges made of stycast epoxy with various angles 
to mount the samples under angled magnetic field in mag¬ 
netization measurements. For the Hall measurement, we 
checked the uniformity of the current flow by using mul¬ 
tiple sets of electric contacts within a sample. The Hall 
data is antisymmetrized in order to remove magneoresis- 
tance component caused by a slight misalignment of the 
contacts. 

The measurement is configured such that the flow of 
current lies within the crystallographic a 6 -plane, while 
the saturated magnetization sweeps the out-of-plane po¬ 
lar angle from 6h = 0 to 90° (see the inset of Figl2] (a)). 
This corresponds to the out-of-plane anisotropic magne¬ 
toresistance configuration as shown in [^,[ 23 . It is dis¬ 
tinct from the typical AMR configuration where both 
current and the magnetization lie within the film plane 
0 . Note that the direction of current always remains 
within the a 6 -plane and therefore, when Oh = 90° the 
current is parallel to H. 

First, we plot the in-plane resistivities (p) as a func¬ 
tion of H for Oh = ^ (thick lines) and Oh = 90° (thin 
lines) at different fixed temperatures (T) in Fig. [TKa) 
and (b). Figlljc) summarizes the T dependence of the 
in-plane magnetoresistance (MR), defined as Ap/po = 
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FIG. 1: (a,b) In-plane resistivity p as a function of H in 
two different orientations: Oh = 0° (thick lines) and Oh = 
90° (thin lines) at fixed Ts (a) T < 50 K and (b) T > 75 K. 
At each T, is clearly visible in both panels but in (b) Hp^ 
is hard to identity because of the large p-axis scale, (c) The 
normalized in-plane MR, Apjpo = [p{H) — Po]/po with po = 
p{H = 0) as a function of T measured at poH = 7 T. Tc = 
133 K is marked with an arrow. Blurred blue line is centered 
at Ta = 52 K, which is determined from the temperature at 
which piPn = 0,77°) is no longer smaller that p(0h = 90°). 
However, since we intend Ta to indicate approximate T scale 
above which the difference between the two MRs diminishes, 
we refer to Ta « 50 K in the main text. 


[p{H,0h) - Po]/po with Po = p[H = 0), at poH = 7 T, 
where all of the spins are polarized along either orienta¬ 
tion of H. The discrepancy of the in-plane MR between 
the two orientations is most striking at 2 K, where the 
reduction of MR with out-of-plane magnetization is three 
times larger than with in-plane magnetization. We will 
come back to this below. Also the 77-independent re¬ 
sistivity values in 77 > 77° and 77 > 77“^ for both 77 
directions imply that all spins are polarized and the con¬ 
tribution of spin-disorder induced scattering is very little. 

As T increases, this trend reverses around Ta « 50 
K, and p{H,0h = 0°) becomes higher in the entire 77 
range. Ta is empirically determined from the temper¬ 
ature at which p{0h = 0 , 77p is no longer smaller than 
p{0h = 90°), as shown in Fig. [He). In T > Ta, the resis¬ 
tivity discrepancy between two spin polarization is over¬ 
whelmed by spin-disorder scattering induced resistance. 
This is also clear from the 77 dependence of two resistiv¬ 
ities are parallel to each other in T > in Fig. [Ijb). 
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The out-of-plane polarization field is clearly iden¬ 
tified up to Tc, as is although the large x-axis scale 
makes it hard to identify in Fig.jljb). At any given T, the 
lower polarization field for 9h = 90° indicates a more ef¬ 
fective suppression of spin-disorder for the same strength 
of H. This also gives rise to a lower p value in entire H 
range for high Ts. The similar rate of reduction of p at 
higher T, which makes the two curves at each tempera¬ 
ture almost parallel, implies that their H dependence is 
attributed to the suppression of spin disorder. 



FIG. 2: (a) The magnetization component parallel to the ap¬ 
plied magnetic field (M||) as a function of H at different angle 
9hs a.t T = 5 K. The measurement configuration is shown in 
the inset, (b) In-plane MR as a function of H oriented at 
angle 9h and T = 5 K. Vertical lines mark the polarization 
field (Hp^) at a given 9h- At intermediate angles, all spin are 
jxfia rized at Hp^ but are not yet aligned along H direction 


In Fig[2] (a) we plot the magnetization along the H 
direction (M||) as a function of H for different ^jjs, at 
T = 5 K <C Tc- The polarization fields (Hp^) are clearly 
visible, indicated with vertical lines and with arrows for 
both end angles oi 9 h = 0 and 90°. It occurs at fields 

well approximated by ~ rapid rise 

in magnetization observed at large angles e.g. 9h = 90° 
or 27° indicates a phase transition from the chiral soliton 
lattice to the ferromagnetic ordering. At intermediate 
angles, the rise at Hp^ is not as rapid as at 9h = 90°. 
While all spins are polarized when H = Hp^ at these 
angles, they are not yet aligned to the direction of H, 
but remain closer to the a&-plane. Upon increasing H 
more, the polarized spins collectively rotate until they 


eventually align to H. This results in the gradual rise 
of M|| in 77 > Hp^ toward the saturated value, which 
is ascribed to the competition between Zeeman and the 
magnetic anisotropy energy [25 1. 

In Fig. 121 (b), we compare this with the in-plane MR. 
At 9h = 90°, the reduction in MR by 5.6% indicates 
a decrease in spin scattering as the helical ordering be¬ 
comes polarized along 77. A clear kink denotes the spin- 
polarization field Hp^. This dome-like shape is com- 
monly observed in other helimagnets for the same reason 
2^ 271. Once the spins completely polarize within the 


plane (77 > 77“°), there is little variation in MR with 77, 
which ensures that the contribution of the spin-disorder 
induced MR is negligible. It is worth pointing out that 
Ap/po at 9h = 90° bears great similarity, in both 77 de¬ 
pendence and magnitude, to that of the interlayer MR, 
i.e. the MR when 7 is applied along the c axis and 77 is 
in the ab plane, reported in 211. 

When 9h approaches 0, upon increasing 77, the MR 
decreases much slower initially but eventually surpasses 
the in-plane value, reaching down to 14.1%, almost 3 
times more reduction than 9h = 90°. This reduction is 
remarkable compared to the variation found in a typical 
traditional AMR phenomena, which are only a few tenths 
to a couple of percent at most [i^. Although slightly 
larger than a couple of percent out-of-plane AMR values 
have been reported in magnetic thin films systems [j^ . 
it is attributed to the reduction of spin scattering due to 
the geometrical size and texture of films. It is only very 
recent that the interface effect between adjacent magnetic 
and non-magnetic film was identified for a possible source 
of the out-of-plane AMR effect, yet the reported size is 
less than 1% 

For 9h = 8° and 12° there is another kink in the MR, 
appearing in the range of 77 where the phase transition 
from chiral soliton lattice to ferromagnetic ordering oc¬ 
curs [2l|. This finite width of the transition, only oc¬ 


curring for angles close to the c-axis, is believed to be 
caused by the formation of multiple domains undergoing 
the transition at slightly different fields, as also seen in 
M||(77) of Fig|2| (a). While these kinks indicate a great 
sensitivity of the electrical transport to changes in the 
spin structure, it does not account for the continued de¬ 
crease in MR in 77 > 77®^. 

To have a better understanding of this spin-orientation 
dependent in-plane MR, we examine the AMR effect in 
both in-plane and out-of-plane rotation of applied field. 
Schematics of the measurement are shown in insets of 
FigEl In the main panel, both in-plane and out-of-plane 
angular dependence of the in-plane resistivity normalized 
with p|| = p{pM = 0, pqH = 7 T) are plotted, where pM 
denotes the angle between M and 7 (see Supplementary 
Material). In both cases, rotation of the magnetization 
give rise to a cos^ pM angular dependences indicated with 
solid lines. However, the normalized AMR oscillation 
amplitude for the out-of-plane MR is larger by 23 times 
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FIG. 3: In-plane MR at = 7 T, normalized with 

P|| (<Pm = 0, fioH — 7 T),where ipM denotes the angle between 
the current and M both in-plane (open) and out-of-plane 
(closed circle), as shown in the sketches. Note the amplitude 
of out-of-plane AMR is 23 times larger than in-plane AMR. 
The angular dependences for both cases are found oc cos^ tpM, 
of which is displayed in the solid lines. Inset is plotted in-plane 
AMR with magnified y-axis. 


than the in-plane one; -1-0.4% and -9.1% change relative 
to p|| for in-plane and out-of-plane respectively. Note 
the different sign indicates the most reduced resistivity 
occurs when M || / for in-plane and M _L / for out- 
of-plane. This out-of-plane AMR effect should also be 
distinguished from giant magnetoresistance of magnetic 
multilayers, in which the saturated MR remain the same 
for both orientations (^ . 

Finally, we examine the Hall effect which is measured 
with H applied along c-axis {9h = 0). This configuration 
matches those in which a topological Hall effect (THE) 
was observed in other helimagnets, e.g. with H along 
the helical axis [l^. Fig. 0](a) shows field sweeps of the 
Hall resistivity {pyx) at various fixed Ts. At low tem¬ 
peratures, an unexpected H dependence of pyx emerges, 
characterized by the pronounced concave curvature for 
0 < H < Hp. As shown in FigHJ^c), this peculiar de¬ 
pendence is more obvious after subtracting the ordinary 
Hall contribution, as described in the forthcoming anal¬ 
ysis. Note that this behavior emerges below Ta « 50 
K, coinciding with the emergence of the MR discrepancy 
described above. 

For T > Ta, the Hall signal exhibits archetypi¬ 
cal anomalous Hall behavior, empirically expressed as 
Pyx{H) = pqRhH -I- pqRsM, where /tq is the vacuum 
permeability, Rh normal Hall coefficient, and Rs the 
anomalous Hall coefficient. The first and second terms 
correspond to the ordinary Hall effect (OHE) and the 
anomalous Hall effect (AHE), respectively. As T de¬ 
creases, the AHE signal becomes visible at T = 150 
K, and grows more prominent through Tc, where it is 
marked by a sharp kink at H = //“. Erom the shape 
of the curves Rs has the opposite sign as Rh- In the 
limit of the intrinsic AHE, we can rewrite the anomalous 
Hall coefficient as Rg = ShP^{T,H), with Sh an H in- 



FIG. 4: (a) H dependence of pyx at various Ts below (solid 
lines) and above (broken lines) Tc- Peculiar concave H- 
dependence occurs only in T < Ta ~ 50 K. Since H || c, 
all spins get smoothly polarized along c-axis via helical and 
conical state as H increases, (b) The ordinary {Rh) and 
anomalous (Sh) Hall coefficients, determined from fits of pyx- 
Above Ta, both Rh and Sh are independent of T. Below Ta, 
the ht no longer works and Sh cannot be obtained in a re¬ 
liable way. Rh is extracted from poH > 3 T where pyx is 
linear to H. (c) H dependence of pyx less the OHE. AHE 
contribution decreases as T is lowered and is replaced with 
the hump-like feature in T < Ta, which has very little re¬ 
semblance with M{H). The curve for 50 K displays a little 
dip before becoming flat, displaying the crossover between the 
hump-like shape and AHE component. 


dependent constant [2^, and fit our pyx{H) traces well 
in Ta < T < Tc- The fitting parameters Rh and Sh are 
plotted in Fig[l](b). The effective carrier density derived 
from Rh in T > Ta is found 9 x 10^° holes/cm^. It is 
consistent with the value found above Tc at 200 K, where 
the Hall signal recovers almost a linear H dependence. 


In T < Ta however, this analysis is no longer valid and 
cannot capture the concave H dependence, even when 
an additional extrinsic contribution term is included. As 
seen in FigH] (c), the emergence of a novel H dependence 
below Hp replaces the gradual reduction of the AHE con¬ 
tribution. It is hard to connect it with the AHE picture, 
as it bears no resemblance to M{H). We rule out the 
possibility that this behavior is caused by the topologi¬ 
cal Hall effect (THE) for the following reasons: first 
our spin structure model found the skyrmion density is 
zero in this material for any temperature and the forma¬ 
tion of skyrmions is unlikely due to the large magnetic 
anisotopy [1^. Second, the T dependence - occurring 
at low Ts and vanishing as T increases - is the oppo¬ 
site to what was observed in the THE in other magnetic 
systems [27|. Lastly, a recent neutron scattering study 
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could not identify the existence of a complex periodic 
spin structure of skyrmions (^ . 

It is interesting to compare the Hall effect of 
Cri/ 3 NbS 2 to that of Fei/ 4 TaS 2 , an anisotropic mag¬ 
netic dichalcogenide with a similar crystalline structure. 
Both have the same resistivity at 5 K, but Cr 4 / 3 NbS 2 has 
one-tenth the carrier density of Fei/ 4 TaS 2 . This means 
the Cr 4 / 3 NbS 2 has a ten-fold larger mean free path (1). 
The intrinsic AHE signal scales with 1/P and diminishes 
rapidly as T is lowered and p decreases [2^, resulting 
in the OHE dominating the Hall signal at low temper¬ 
atures, i.e. a recovery of H linear dependence of pyx- 
This is consistent with our observation that the extrinsic 
AHE is ne glig ible over the entire temperature range in 
Cri/3NbS2M. 

The fact that unusual H profile of pyx and the larger 
amplitude of out-of-plane AMR occur in the same T 
range, T < Ri 50 K, implies that the two phenom¬ 
ena share a similar origin. In this temperature regime 
Rh, which was estimated from the slope of pyx when 
H > 3 T, increases rapidly. The crossover temperature 
Ta also appeared in other measurements: the T depen¬ 
dence of the thermopower (see Supplementary Material) 
changes slowly just below Tc, and then after a broad 
shoulder around 50 K rapidly decreases. These unusual 
spin-orientation dependent transport features, unique to 
Cri/ 3 NbS 2 , emerge only below Ta- Ta is also consis¬ 
tent with a deviation from Bloch’s T^P law of M(T)[^. 
These observations points to a spin-orbit coupling effect 
of order ksTA, with ks the Boltzmann constant. 



FIG. 5: The calculated magnetic density-of-states of 
Cri/ 3 NbS 2 for the indicated magnetic moment orientations; 
(a) total DOS, (b) Cr DOS. The energy scale ± 25 meV cen¬ 
tered at Fermi energy Ep is magnified in the insets. The 
DOS at Ef was found to be 3.2 % greater when the moments 
orient along the c-axis instead of perpendicular to it, while 
the in-plane plasma frequencies are essentially unaffected by 
magnetization orientation, which is consistent with enhanced 
conductivity within an order of unity. 



E-Ep (eV) 

FIG. 6: The calculated magnetic density-of-states of 
Gri/ 3 NbS 2 projected by spin and atomic character. Inset: 
the orbital-projected character of the Nb DOS. 


In an attempt to understand these anisotropic magne¬ 
totransport properties, first principles calculations were 
pursued (procedural details in the Supplementary Mate¬ 
rial). Results of these calculations, depicted in Fig. 5 
show a significant dependence of the calculated magnetic 
density of states on moment orientation, both for the to¬ 
tal DOS (top) and Cr DOS (bottom). As indicated in the 
inset to Fig. 5, for states very near the Fermi level (Ej?) 
- less than 25 meV from E^t- - the DOS can vary by as 
much as 15% depending on moment orientation. More¬ 
over, the Fermi level density of state (DOS) was found 
to be 3.2 % greater when the moments orient along the 
c-axis instead of perpendicular to it, while the in-plane 
plasma frequencies are essentially unaffected by magne¬ 
tization orientation. Thus, neglecting spin-orientation 
dependent scattering processes, the increase in DOS is 
largely responsible for the reduction of the in-plane resis¬ 
tivity when the magnetization lies parallel to the c-axis 
instead of the a&-plane. The resistivities for the two spin- 
polarization directions at 5 K differ by 11%, which is big¬ 
ger than the DOS difference by a factor of the order of 
unity. Also, by examining the bands near Ep (within a 
20 meV window) for both orientations, one finds the aver¬ 
aged offset in the energies to be approximately 1.7 meV 
cs 20 K. This is a direct result of SOC and consistent 
with the temperature scale Ta found in the MR and Hall 
effect data. These changes in electronic structure and 
their subsequent modification of scattering processes are 
likely responsible for the spin orientation dependence in 
the transport. 

We note also from Fig. 5 that the Cr-site DOS is of 
essentially the same shape as the total DOS, within an 
eV of Ef. In addition, the effects on the Cr DOS of the 
change in moment orientation from [001] (|| c) to [100] (|| 
ab) closely parallel the changes in the total DOS. This is 
due to the primary role of the Cr atoms in the magnetic 
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behavior - only for atoms with a large local moment, such 
as the Cr here, does one expect a moment-orientation 
dependent DOS. Note that in the calculation the induced 
Nb moment is only 0.05 and that on S essentially zero. 

In Figure 6 we show the magnetic state DOS projected 
into spin-up and spin-down sub-bands, along with the Nb 
Ad and 5s orbital-projected DOS. A substantial exchange 
splitting is visible in the Cr DOS, as expected given its 
magnetic character, with some minor splitting also ob¬ 
served in the Nb and S. There is little, if any, spin-down 
character to the Cr DOS for nearly 2 eV below Ef. This 
is suggestive of the strength of the magnetism in this 
system, as indicated by the 3 /is/Cr total magnetic mo¬ 
ment found both in the calculation and the experiment. 
There is substantial hybridization around indicative 
of covalent bonding, although the Cr states are generally 
confined to within 2 eV of Ef and the Nb and S predom¬ 
inate outside this range. In the inset of Figure 6 we plot 
the orbital-projected Nb DOS, which is dominated by Ad 
contributions, with the 5s contribution negligible. 

Further detailed investigation on the electronic struc¬ 
ture depending on the spin direction will shed light on 
understanding the difference of amplitudes of in-plane 
and out-of-plane AMR. This mechanism is suspected to 
be closely related to abnormal H dependence of the Hall 
effect observed for T <Ta- 

The remarkable sensitivity of the magnetotransport 
behavior to polarized spin orientation in Cri/ 3 NbS 2 is 
mostly driven by Cr^^ ions. Strong magnetization orien¬ 
tation dependence of the electronic structure of Sd met¬ 
als leads to large effects on the magnetotransport prop¬ 
erties in Cri/ 3 NbS 2 . The same strong dependence also 
underlies efforts to control magnetic anisotropy in mul¬ 
tilayer systems made of id magnetic elements. Our re¬ 
sults demonstrates a closely interwoven relation among 
the orientation of magnetization, the resultant electronic 
structure change and its influence on the electrical con¬ 
duction, all of which are mediated by spin-orbit coupling. 
Furthermore, the non-collinear helical ground state en¬ 
riches the role of spin-orbit coupling and the subsequent 
complexity of the electrical transport properties. 

Acknowledgment The authors thank J. G. Checkelsky, 
X. Fan, M. Hermele and K. M. McElroy for enlightening 
discussions. This work was supported by the US DOE, 
Basic Energy Sciences, Materials Sciences and Engineer¬ 
ing Division (ORNL) and at CU under Award Number 
de-sc0006888, D.P. is supported by the Critical Materi¬ 
als Institute, an Energy Innovation Hub funded by the 
U.S. DOE, Energy Efficiency and Renewable Energy, Ad¬ 
vanced Manufacturing Office. 

Current Address: Los Alamos National Laboratory, Los 
Alamos, NM 


* Electronic address: minhyea.lee@colorado.edu 

[1] H. Brooks, Phys. Rev. 58, 909 (1940), URL 

http://link.aps.org/doi/10.1103/PhysRev.58.909 

[2] P. Bruno, Phys. Rev. B 39, 865 (1989), URL 

http://link.aps.org/doi/10.1103/PhysRevB.39.865 

[3] D.-s. Wang, R. Wu, and A. J. Free¬ 
man, Phys. Rev. B 47, 14932 (1993), URL 

http: //link, aps . org/doi/10.1103/PhysRevB. 47.14932 

[4] M. Bode, S. Heinze, A. Kubetzka, O. Pietzsch, 

X. Nie, G. Bihhnayer, S. Bliigel, and R. Wiesen- 
danger, Phys. Rev. Lett. 89, 237205 (2002), URL 
http://link.aps.org/doi/10.1103/PhysRevLett.89.237205 

[5] F. Gimbert and L. Galmels, Phys. 

Rev. B 86, 184407 (2012), URL 

http://link.aps.org/doi/10.1103/PhysRevB.86.184407 

[6] K. Hotta, K. Nakamura, T. Akiyama, 

T. Ito, T. Oguchi, and A. J. Freeman, 

Phys. Rev. Lett. 110, 267206 (2013), URL 

http://link.aps.org/doi/10.1103/PhysRevLett.110.267206 

[7] S. Muhlbauer, B. Binz, F. Jonietz, G. Pfleiderer, 

A. Rosch, A. Neubauer, R. Georgii, and P. Boni, Science 
323, 915 (2009). 

[8] M. Lee, W. Kang, Y. Onose, Y. Tokura, and N. P. Ong, 

Phys. Rev. Lett. 102, 186601 (2009). 

[9] A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, 

R. Ritz, P. G. Nikfowitz, and P. Boni, 

Phys. Rev. Lett. 102, 186602 (2009), URL 

http: //link, aps . org/doi/10.1103/PhysRevLett. 102.186602 

[10] I. Dzyaloshinsky, Journal of Physics and Ghem- 

istry of Solids 4, 241 (1958), ISSN 0022-3697, URL 

http://www.sciencedirect.com/science/article/pii/0022369758 

[11] T. Moriya, Phys. Rev. 120, 91 (1960), URL 

http://link.aps.org/doi/10.1103/PhysRev.120.91 

[12] J. M. Hopkinson and H.-Y. Kee, Phys. 

Rev. B 79, 014421 (2009), URL 

http://link.aps.org/doi/10.1103/PhysRevB.79.014421 

[13] A. Fert, V. Gros, and J. Sampaio, Nature nanotechnology 
8, 152 (2013). 

[14] N. Romming, C. Hanneken, M. Menzel, J. E. Bickel, 

B. Wolter, K. von Bergmann, A. Kubetzka, and 
R. Wiesendanger, Science 341, 636 (2013). 

[15] B. G. Sales, R. Jin, and D. Mandrus, 

Phys. Rev. B 77, 024409 (2008), URL 

http://link.aps.org/doi/10.1103/PhysRevB.77.024409 

[16] Y. Mokrousov, H. Zhang, F. Freimuth, B. Zimmermann, 

N. H. Long, 1. Weischenberg, J.and Souza, P. Mavropou- 
los, and S. Blugel, J. Phys.: Gond. Matt. 25, 163201 
(2013). 

[17] N. J. Ghimire, M. A. McGuire, D. S. Parker, 

B. Sipos, S. Tang, J.-Q. Yan, B. C. Sales, and 
D. Mandrus, Phys. Rev. B 87, 104403 (2013), URL 
http://link.aps.org/doi/10.1103/PhysRevB.87.104403 

[18] T. Moriya and T. Miyadai, Solid State Comm. 42, 209 
(1982). 

[19] T. Miyadai, K. Kikuchi, H. .Kondo, S. Sakka, M. Arai, 
and Y. Ishikawa, J. Phys. Soc. Jpn. 52, 1394 (1983). 

[20] Y. Togawa, T. Koyama, K. Takayanagi, S. Mori, 

Y. Kousaka, J. Akimitsu, S. Nishihara, K. In- 
oue, A. S. Ovchinnikov, and J. Kishine, 

Phys. Rev. Lett. 108, 107202 (2012), URL 

http: //link, aps . org/doi/10.1103/PhysRevLett. 108.107202 



7 


[21] Y. Togawa, Y. Kousaka, S. Nisliiliara, K. Inoue, 

J. Akimitsu, A. S. Ovchinnikov, and J. Kisliine, 

Phys. Rev. Lett. Ill, 197204 (2013), URL [27] 

http://link.aps.org/doi/10.1103/PhysRevLett.Ill.197204 

[22] T. G. S. M. Rijks, S. K. J. Lenczowski, R. Coehoorn, and [28] 
W. J. M. de Jonge, Phys. Rev. B 56, 362 (1997), URL 

http: //link, aps . org/doi/10.1103/PhysRevB. 56.362 

[23] A. Kobs, S. Hefie, W. Kreuzpaintner, G. Winkler, 

D. Lott, P. Weinberger, A. Schreyer, and H. P. [29] 

Oepen, Phys. Rev. Lett. 106, 217207 (2011), URL 

http: //link, aps . org/doi/10.1103/PhysRevLett. 106.2172(p0Ci] 

[24] T. McGuire and R. Potter, IEEE Trans. Magn. 11, 

101838 (1975). 

[25] B. J. Chapman, A. Bornstein, N. J. Ghimire, D. Man- 
drus, and M. Lee, Applied Physics Letters 105, 072405 
(2014). 

[26] M. Lee, Y. Onose, Y. Tokura, and N. P. 


Ong, Phys. Rev. B 75, 172403 (2007), URL 

http: //link, aps . org/doi/10.1103/PhysRevB. 75.172403 
B. J. Chapman, M. G. Grossnickle, T. Wolf, and M. Lee, 
Physical Review B 88, 214406 (2013). 

M. N. Baibich, J. M. Broto, A. Fert, F. N. Van Dan, 

F. Petroff, P. Etienne, G. Creuzet, A. Friederich, and 
J. Chazelas, Phys. Rev. Lett. 61, 2472 (1988), URL 
http://link.aps.org/doi/10.1103/PhysRevLett.61.2472 

N. J. Ghimire, Ph.D. thesis. University of Tennessee 
(2013). 

The second term py^ = pqRhH + poRsM was modified 
to pLoSHp^M {\+ — ), with pi an additional fitting param¬ 
eter, but fits revealed ^ <C 1 over most of the T range 
investigated. Meanwhile, in Fei/ 4 TaS 2 , Pxy is dictated by 
the intrinsic AHE at low T and then the extrinsic AHE 
as T approaches Tc- 


